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We perform some systematic numerical search for Schwinger boson mean field states on the square, 
triangular and kagome clusters. We look for possible inhomogeneous ground states as well as low- 
energy excited saddle points. The spectrum of the Hessian is also computed for each solution. On 
the square lattice we find gapless (7(1) gauge modes in the non-magnetic phase. In the Z2 liquid 
phase of the triangular lattice we identify the topological degeneracy as well as vison states. The 
energy landscape on kagome lattices shows tiny energy scales and we discuss its relation to classical 
planar states. 



I. INTRODUCTION 

Quantum antifferomagnets can display complex many 
body phenomena, with rich phase diagrams, exotic states 
of matter with emerging degrees of freedom. Indeed, 
minimizing a rather simple looking interaction like the 
Heisenberg one, Si ■ Sj, can lead to a vast variety of 
states of matter, depending on the size Sf = S(S + 1) of 
the spins, the space dimension, geometry of the lattice, 
or relative strengths of possible competing interactions 
(frustration, etc.). In fact, after many years of theoreti- 
cal investigations, the nature of the ground state of the 
spin- 1 Heisenberg model remains controversial on several 
two and three dimensional lattices (like the kagome and 
pyrochlore) . Some of the most interesting states that can 
be stabilized at T — are called spin liquids and have 
no direct classical analogs. These systems remain rota- 
tionally invariant down to zero temperature. There is a 
whole zoo of possible spin liquids, and the most exotic 
one have low energy excitations which carry a half-odd- 
integer spin (unlike conventional spin waves). 

There are rather few theoretical tools that are able to 
describe the limit of strong quantum fluctuations (small 
S) for these systems, and spin liquid phases in particular. 
Among them, the so-called large-iV techniques, intro- 
duced by Affleck for spin chains, J play a central role. By 
generalizing the symmetry group of global spin rotations 
from SU(2) to some larger group like SU(jV) or Sp(27V), 
the model can be solved in the limit N = 00. Of course, 
the physics for N = 2 needs not be simply related to that 
at N = 00. However, the success of these approaches is 
in part due to the fact that a number of interesting states 
that can be realized for SU(2) models, also exist in the 
N = 00 phase diagrams. In particular, the N = 00 
models are not restricted to have magnetically ordered 
ground state, and can have resonating valence-bond spin 
liquids ground states with fractionalized excitations and 
emerging gauge degrees of freedom. Furthermore, the 
effect of Unite-N corrections (both perturbative 4 ' ' and 
non-perturbative ones ) can be addressed. 

In this work we consider a particular large-A^ 
limit, the so-called Schwinger boson mean field the- 
ory (SBMFT). 1 ' 1 It has been applied to many two- 



dimensional systems, such as the J\ — J 2 square, 
triangular, ~ honeycomb 12 ' 1 ', kagome, i,,14_1<, ) kagome 
with further neighbor interactions, ' or Shastry- 
Sutherland 1 v ' or CaV4C>9, ,! lattices. This limit can 
in particular stabilize magnetically ordered (Neel) states 
as well as Z 2 spin liquids with gapped spinous." 

The model contains N flavors of spin-i bosons 
(spinons) and the parameter which plays the role of the 
spin value 25* is the number k of bosons per site (and per 
flavor) . After performing a Hubbard-Stratonovich decou- 
pling of the boson-boson interaction, the different boson 
flavors are coupled to a single complex field Aij on each 
bond ij of the lattice. The formal (Gaussian) integration 
of the spinons gives an effective action iSeff.jv^A] for the 
bond field. But the dependence in N is simple since this 
action is that of the single-flavor problem multiplied by 
a factor N: S cS , N [A] = NS eS< x[A]. When N is large it 
is therefore natural to perform a saddle point expansion. 
At the saddle point the bond- field fluctuations are frozen, 
and the effective Hamiltonian is simply quadratic in the 
spinon operators. So the action can be computed with 
the help of a standard Bogoliubov transformation. Find- 
ing a the large- N ground state thus amounts to solving 
a classical minimization problem. 

Previous SBMFT studies have mostly been focused on 
the ground state properties at N = 00. However, as 
a first step to understand the effects of finite- N correc- 
tions, it is interesting to have access to the low-energy 
energy landscape of saddle points. Although the system 
is locked into the lowest one when N — 00, the (large 
but) finite- TV physics must include some fluctuations be- 
tween different low-energy saddle points, as well as some 
(perturbative in 1/N) fluctuations in the vicinity of each 
of these saddle point. Our goal here is to provide a quan- 
titative description of the/some excited saddle points. 

In order to reduce the number of variables, it is almost 
always assumed that the lowest-energy state preserves 
most (or at least some) of the lattice symmetries (see 
Ref. 21 for a notable exception). In practice most numeri- 
cal studied have so far been restricted to a solutions with 
a unit cell including a few sites only. This is a reasonable 
assumption for the ground state, and we confirm valid- 
ity in several cases. However, since we are interested 
here in excited saddle points, we need to be able to com- 
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pute some spatially inhomogeneous states as well. We will 
however limit ourselves to time-independent solutions. 

In this work perform some extensive numerical min- 
imization to look for the ground state and low-energy 
excited states on clusters containing up to 144 sites. We 
compare the results we obtain on three different lattices 
(square, triangular and kagome). and for different values 
of the boson density ("spin") k. 



II. SCHWINGER BOSON MEAN FIELD 
THEORY 

To keep this article self-contained, this section presents 
the basic ideas and notations of the SBMFT (see also 
Refs. 22 ' 23 ). 

"up" and "down" bosons operators (a € {T, 4-}) j carry- 
ing a S = 1/2, are introduced at each lattice site: b\ a 
and b ia . The spin operators can then be written: 



?+ _ uf 



bLb. 



25* = b\ t b it -blb u 



(1) 
(2) 
(3) 



These relations imply that the commutation relations 
[5f ,Sf] = ie a/3S Sf are automatically verified. The total 
spin reads Sf — ^ + l), where m = b\^.b^ + b\^bn is 
the total number of bosons at site i. To fix the "length" 
of the spins, the following constraint must therefore be 
imposed on physical states: 



V», Y>blb Ur = K = 2S 



(4) 



The Heisenberg exchange Hamiltonian is biquadratic 
in the b operators and reads 



H = ^\/„S. -S, (5) 



where a is the vector whose components are the Pauli 
matrices, and each lattice bond ij is taken only once. It 
is convenient to re write H using rotationally invariant 
bond operators: 



(7) 



A\- creates a spin singlet on the (oriented) bond ij 
whereas Bl- creates a triplet. Due to the constraints, 



these operators are linked by the relation 

5 2 . 



bIb, 



-At; A 



(10) 



and it is therefore possible to express the Hamiltonian 
using A only: 



H 



5> 

(ij) 



s- 



2A^Aij 



(11) 



The SB mean field approximation - which can be for- 
mally justified in a large- N limit of the model- consists 
in decoupling the quartic terms and to add some chemi- 
cal potentials Aj to tune the average number of boson at 
each site site (instead of Eq. 4) . The resulting mean field 
Hamiltonian is 



Ih-.v ■■ > -A , ■( 5 2 - 2A ij A ij - 2A ij A\ j + 2 



5>( 



(ij) 



(12) 



where x is the conjugate of x, Aij are complex link vari- 
ables with property Ay — —Ay L . It is also possible to 
write a mean field theory keep simultaneously both oper- 
ators A and B on each bond. 24 ' 25 This will however not 
be considered here. 26 

It is convenient to write Bmf by grouping the cre- 
ation and annihilation operators in a vector <f> = 



(b n ,---,b Nt ,b\ v ---,tf Ni y : 

H M f = (4>) f M$ + est. 
where the M is a 2N x 2N matrix : 



M = 



J ij A-ij 



JijAij 

-A,; 



(13) 



(14) 



For simplicity, we restrict the discussion to finite systems 
where Hmf nas a g a PP e d spectrum (possibly vanishing 
in the thermodynamic limit). The ground state exists 
and its spectrum is gapped if and only if the the matrix 
M is positive definite. We define a diagonal 2N x 2N 

q ^ and, with the gap condition, 

crM is diagonalizable (although not Hcrmitian) and has 
pairs of real eigenvalues ±w„=i - at, where the u> n > are 
the of the Bogoliubov modes of Hmf ■ The ground state 
energy of Hmf is given by 



matrix a: a = 



(ij) 



E 



MF — 



= £ w » + Yl J v i^j i 2 + ( 25 + 1) E A « ( 15 ) 



where : — : represents normal ordering and the bond 
operators Ay and -By are 



A 



(8) 
(9) 



(remark: in a large- N formalism, this corresponds to the 
energy per flavor) . 

The self consistency is reached when Ay =< A^ >gs, 
which is equivalent to dE/dAij — 0, where E is the en- 
ergy of the ground state of Hmf- The average num- 
ber of boson per site is equal to 25 (Eq. 4) at the 
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point(s) where dE/dXi — 0. Since d 2 E/d\ 2 < (see 
this footnote 2 ^), this point corresponds to a maximum of 
{Xi} H> E[{Aij} , {Aj}]. For fixed this allows to 

use a maximization algorithm to determine the chemical 
potentials Ai[{^4ij}]- In a similar way, the self-consistent 
correspond to a minimum of E [{Aij} , A^j^j}]]- 



A. Gauge invar iance and fluxes 

The solution of this system of equations is not unique, 
at least because of the U(l) gauge invariance of the 
Hamiltonian. Under the gauge transformation 

bja _> e i8 ^b ja 

A tj -> jWWWAij (16) 

the physical operators (as Si or H) stay unchanged. In 
other words, the bond parameters label the physical 
state in a redundant way. Using gauge transformations 
we can fix the phases of some A^ , without changing any 
physical observable. This decreases the number of vari- 
ables to be optimized, and is therefore useful in numerical 
studies. 

The moduli \Aij\ are of course gauge invariant quanti- 
ties, simply related to the energy. But some combinations 
of complex phases around closed loops are also gauge in- 
variant. Consider the following operator Aj 1 j 2 ...j 2n : 

&iii 2 -i 2 n — Ai^zi— A\ 2i:i )A l3U ...(— A i2nii ). (17) 

This operator is defined on any loop with an even 
length on the lattice and is manifestly gauge invariant. 
Its mean field (or large- N) counter part 

gives the flux 4>i 1 i 2 ---i 2n = ar g(-^iii 2 ---i2n) recently dis- 
cussed by Tchernyshyov et al.. 

In general, the number of variables (all the moduli, 
plus some phases) grows with the system size. Since 
it is quite difficult to perform an exhaustive search for 
mean field solutions if the number of parameters is ex- 
tensive, the usual strategy is to decrease the number of 
bond variables by assuming that the ground state solu- 
tion preserves some (or all) symmetries of the lattice. We 
will however give some examples where the lowest energy 
solution spontaneously breaks some symmetries. 



III. NUMERICAL METHOD 

We describe an algorithm to find saddle points (self- 
consistent mean field state), local minima, and (hope- 
fully) global energy minima. 



A. Determination of the chemical potentials 

For given values of the bond parameters Ay, one 
should first adjust the chemical potentials (Aj). 

The basic idea is to perform a (non-linear) least-square 
minimization of where /j({Aj}) = (hi) — k is a 

function of the chemical potentials. Since each density 
(hi) is an increasing function of A^, this method con- 
verges relatively rapidly. To do so, we use an imple- 
mentation of the Levenberg-Marquart algorithm. To 
make the method faster, we provide the gradient ma- 
trix Gij = explicitly. Gij is easy to compute us- 
ing linear-response theory (using the matrix which im- 
plements the Bogolibov transformation, as well as the 
energies Wj)). The iterative minimization should start in 
a region of the space of A where the Bogoliubov trans- 
formation exist (M > 0). If the previous values of A do 
not satisfy this condition with the new one starts 
the least-square minimization of ff from a sufficiently 
low uniform A. 

Finally we note that for some values of the bond pa- 
rameters Aij, there is no uncondensed state satisfying 
(hi) = k. If such a situation is encountered, we add some 
artificial energy penalty so that the energy-minimization 
algorithm (next section) tends to escape this point. 

B. Optimization of the bond parameters Aij 

Equipped with a procedure to compute the A.; as 
a function of the A^. we can start to look values 
of A^ which correspond to self-consistent mean field 
states. The first stage is simply to iterate the bond self- 
consistency conditions in the usual way: 

We start from initial random bond values (or per- 
turbing a previously found solution) . 

1 After adjusting the A by the procedure above, the 
ground state is obtained by Bogoliubov transfor- 
mation, and the expectation values (A^) are com- 
puted. 

2 The bond parameters are replaced by the val- 
ues above. The new bond parameters are gauge- 
transformed to a fixed gauge choice where a maxi- 
mum number of bond parameters are set to be real. 
This avoids some possible slow drift of some com- 
plex phases of the Aij which would be non-physical. 

3 Go back to step 2 until the bond parameters do not 
change by more than a small threshold e. 

This method was used by Hermele et al. in a similar 
context. It is easy to check that local minima are attrac- 
tive for these iterations while local maxima a repulsive. 
More precisely, the error will decrease (resp. increase) in 
the directions corresponding to eigenvectors of the Hes- 
sian (anticipating on Eq. 20) with positive (resp. neg- 
ative) values. The convergence is faster if the Hessian 
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eigenvalues are large and positive. It practice, the itera- 
tions above allows to quickly go down in energy, but it is 
not efficient to achieve a full convergence when the sys- 
tem has more than ~ 10 bonds or so. Indeed, it turns out 
that a high accuracy is required to resolve the possibly 
small differences (bond modulations, etc) between differ- 
ent saddle points (kagome case in particular, Sec. IV C). 

To find accurately self-consistent state, the second 
stage amounts to perform a least-square optimization on 
J2ij9ij, where gij({A}) = (Ay) - Ay is a function of 
the bond parameters. Of course, the expectation values 
(Ay) are computed after the Ai({A}) have been deter- 
mined. This second stage allows to converge not only to 
local extrema but to saddle points as well. Again, we 
use a Levenberg-Marquart minimization algorithm, with 

explicit calculation of the gradient d g^^ ■ The later is 
again obtained using linear-response theory, with the ad- 
ditional complexity that it contains some terms coming 
from the variation of the A : 



d(Ai 



dA k 



kl 



dA 



kl 



A fixed 



d{At 



a . b 



d\ a 



d(h) 



ab 



djfib) 
dA 



ki 



(19) 



As for the first stage, we work with bond variables 
corresponding to a fixed gauge choice. In particular, the 
matrix of Eq. 19 is evaluated in the subspace of bond 
parameter variations which is orthogonal to pure gauge 
transformations. This is important to get rid of unphys- 
ical slow phase drifts during the iterations. Most of the 
time a double precision accuracy in reached in less than 
10 Levenberg-Marquart iterations when the system has 
less than one hundred bonds or so. 

Finally the optimization is repeated using at least a 
few hundred (often thousands) of random initial con- 
ditions for the A^. This is relatively time consuming 
since each evaluation of the energy (or its derivatives) 
requires a numerical adjustment of the chemical poten- 
tials (Sec. Ill A), which is itself a (convex) least-square 
problem with many variables. 

Varying the parameter e (stopping criterion for the it- 
erations of stage 1) allows to tune if the method will 
converge toward very low-energy saddle points (small e), 
or saddle point at higher energy (larger e). If e is very 
small (say 10~ 5 ) the first stage will terminate close to 
the ground state and the second step (least-square) will 
converge to the ground state with high probability if the 
system is not too large (< 100 bonds). On the other 
hand, if e is too large, the first stage will stop at some 
relatively high energy configuration. In such region of 
the A space we expect a very (exponentially) high den- 
sity of saddle points. In practice this density of saddle 
point is so high that the program will find a new solu- 
tion at every run, and a given saddle point will rarely be 
obtained twice. The best choice is to adjust e so that the 
iteration stage leads to configurations in a typical energy 



range above the ground state where the number of sad- 
dle points is not too large (a few tens). In such a case, 
after a sufficiently large number of runs, one obtains the 
full (or almost full) list of saddle points in that energy 
window. Of course, due to the large number of variables, 
one cannot exclude the presence of some additional sad- 
dle points with a small basin of attraction with respect 
to this algorithm. 



C. Hessian and stability 

The least-square procedure described above leads to 
self-consistent mean field states, satisfying Ay = (Ay). 
These states are saddle point of the energy (considered 
as a function of the Ay , and the chemical potential being 
themselves functions of the Ay). To check if each mean 
field state is a local minimum, local maximum, or generic 
saddle points with stable as well as unstable directions, 
we compute the Hessian matrix: 



d 2 E 



MF 



1.1' 



dAfdAf, 



4Ji 5 w 5 e . 



d(Af 



(20) 



(n) —k i 



where I and I represent two bonds and e and e denote the 
real or Imaginary part of the bond variables (see Eq. 19 
for derivatives of the bond expectation values at fixed 
densities). 

Due to the gauge invariance of the model, the Hessian 
contains some zero eigenvalues associated to infinitesimal 
gauge transformations. The number of such gauge modes 
can be computed on each given lattice (using the rank of 
a modified adjacency matrix of the lattice ), and these 
non-physical zero eigenvalues are of course be discarded 
when discussing these solutions. For the clusters studied 
here, the number of zero eigenvalues is always equal to 
the number of pure gauge modes, so we can conclude that 
there is no physical zero mode in the Hessian. The sign 
of the smallest non-zero eigenvalue of K tells us whether 
the mean field state is a local minimum, or an unstable 
saddle point. 

The spectrum of the Hessian for the ground state gives 
some information about the magnitude of the 1 /N correc- 
tions due to Gaussian fluctuations in the vicinity of the 
energy minimum. Finally, the spatial structure of the 
lowest eigenvector of the Hessian gives some information 
about the nature of these fluctuations. 



IV. NUMERICAL RESULTS 

A. Square lattice 

The SBMFT phase diagram is well known on the 
square lattice: magnetic long-range order for K > 0.39. 
and a disordered phase with gapped spinons for k < 0.39 
("Coulomb" phase, unstable at finite- AT 1 '). In the ordered 
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phase the spinon gap drops as ~ 1/N S (N s is the number 
of sites) and Bose-condensation occurs in the thermody- 
namic limit, leading to spontaneous break down of the 
spin rotation symmetry. By considering here only finite 
clusters, the ground state is always rotationally invariant 
and the gap finite. Still, both phases can be distinguished 
using standard finite size-scaling for the gap or spin-spin 
correlations. We will discuss how the "energy landscape" 
of mean field saddle points differ between the two phases. 

1. Hessian of the ground state and gauge modes 

We did some extensive search for saddle points on 
square lattices with 36 sites and with k = 2S = 0.1 
and k = 1. The ground state, as expected, is spatially 
uniform, real, and has a vanishing flux on all the square 
plaquettes, whatever the boson density. However, the 
spectrum of the Hessian is quite different in the magnetic 
phase and in the disordered phase. 

In Fig. 1 the smallest eigenvalue of the Hessian (Eq. 20) 
is plotted as a function of the "spin" S, for different sys- 
tem sizes (up to 144 sites). Although finite-size effects 
are important, these data indicate that the Hessian is 
gapped in the thermodynamic limit for large S, while it 
becomes gapless for small S. The transition very likely 
coincides with that of magnetic long-range order. 

The gaplessness of the Hessian in the disordered region 
is due to the bipartite character of the lattice. On a 
bipartite lattice, the bond parameters are invariant under 
staggered gauge transformations: 

b ra — ► br*et- ire (21) 

which, in Wen's terminology, means that the invariant 
gauge group (IGG) is U(l). 

If we perform a spatially varying gauge transformation 
including the staggered factor we get 

b ra — > brafl-W (23) 

Arr , _> ^eWM-'CO) , (24) 

where r (resp. r') is on the even (resp. odd) sublattice. 
Starting from a state described by A® r , we construct a 
phase fluctuation of the bond parameters parametrized 
by a rr > G R: 

A rr j = A rr t -\- dA rT ' == A rr /e rr (^^) 

From Eq. 24, a rr ' transforms conventional U(l) 

gauge field: 

a rr , — > a rr > + 0{r) - 6{r'). (26) 

So, if we denote by E({a rr /}) the energy of the perturbed 
mean field state (after the appropriate adjustment of the 
chemical potentials), the energy should be gauge invari- 
ant under Eq. 26. 



In the small- 5 phase where the bosons are gapped and 
their correlation length is short, E({a rr '}) should be a 
local function of the fluctuation a rr ' . 

The simplest local term compatible with gauge invari- 
ance would be the lattice version of the magnetic energy 
(V x a) 2 , and would take the form of the square of the 
circulation of a around small loops (this can also be ob- 
tained from a small-5 expansion ). This would give an 
Hessian eigenvalue scaling as the square of the smallest 
available wave-vector, that is ~ 1/L 2 ~ 1/N S (L the lin- 
ear size and N s the number of sites). Such modes are 
indeed found in the spectrum of the Hessian and corre- 
spond to the second non-zero eigenvalue. The associated 
eigenvector is displayed in Fig. 2. The thickness of each 
bond I is related to the modulus \dAi \ while the color rep- 
resents the complex argument of dAi/Ai. In the present 
case, these complex arguments take only two values: 
±7r/2 (blue and yellow), indicating that dA is a gauge 
mode of the form of Eq. 25, with dAi/A® ~ iai e iM.. 
In the present case a/ ~ cos(2irxi/L) for a vertical link 
I (oriented from the even to the odd sublattice) at hori- 
zontal position x, and a/ = on horizontal bonds. This 
gauge model may be viewed as a low-energy "photon" of 
the effective gauge theory. 1 ' 

In fact the numerical data indicate that the low- 
est eigenvalue of the Hessian decays faster than 1/N S 
(Fig. 1). The associated mode is represented in Fig. 3, 
and is also a gauge mode (dAi/A® is purely imaginary). 
Inspecting the sign of a/ one sees that it corresponds to a 
change in the "global" flux associated to the large loops 
encircling the torus while local loops are unaffected by 
this gauge mode. The structure of this eigenvector of the 
Hessian turns out to be the same for all system size, and 
does not depend on the spin value. We expect its eigen- 
value to decay exponentially with the system size in the 
small- S phase. 

In the Neel phase, E({a rr '}) need not be short-ranged 
and the arguments above fail. Indeed, the spinon are 
charged particles for the gauge field and their condensa- 
tion gaps out the gauge degrees of freedom (Anderson- 
Higgs mechanism). In the Neel phase, computing the 
boson energy in presence of a global flux (through loops 
encircling the torus) amounts to impose some twist on 
the spin directions and the finite spin stiffness p natu- 
rally leads to an energy cost proportional to the square 
of the flux, and leads to a finite Hessian eigenvalue pro- 
portional to p. In turn, the finite Hessian gap in the Neel 
phase indicates a relative stability of the mean field state 
with respect to Gaussian 1/N corrections. 



2. Excited mean field solutions for k = 1. 

For k = 1 (Neel phase), our search for saddle point 
on the 36 site cluster shows that the ground state 
(—30.93605205126) is well separated from the first ex- 
cited saddle point (E = -28.82425530821, Fig. 4). No 
local minimum was found at low energy (see Tab. I), but 



G 




sq64 E=-6. 7269982 2S=0.2 



FIG. 1: Square lattice model. Top: spin gap as a function of 
k = 2S. In the thermodynamic limit this state is associated to 
magnetic long-range order (vanishing spin gap) for 2S > S c ~ 
0.39. Middle panel: smallest eigenvalue of the Hessian for 
the ground state state on 36, 64 and 144-site square lattices. 
Bottom: scaling of the lowest Hessian eigenvalue as a function 
of the system size N 3 . This shows that the Hessian is gapless 
for 2S = 0.2 and gapped for 2S = 0.8. The Hessian gap is 
presumably finite in the magnetic phase, while it is likely to 
vanish in the thermodynamic limit for S < S c - For 2S = 0.2, 
first eigenvalues decays faster than 1/N a whereas the second 
eigenvalue goes to zero as ~ 1/N 3 (see text). 



only saddle points. In addition, the excited saddle points 
turn out to have very unstable directions (strongly neg- 
ative Hessian eigenvalues). There might be some other 
mean field states in the energy range of Tab. I, but these 
should have rather small basin of attraction with respect 
to our search algorithms, since this list of the twelve low- 
est energy state is stable after thousands of runs starting 
from random initial conditions. 

The first excited saddle point is displayed in Fig. 4. 
It takes the form of an excitation localized around one 
site (here site number 14). This state is chiral, with a 
flux equal to ±0.4046133271 on the four square plaquettes 
touching the center of the excitations. The flux then 



FIG. 2: Hessian eigenvector dAi corresponding to the sec- 
ond smallest eigenvalue (0.26224396339) for a 64-site square 
lattice with S = 0.1. The thickness of each bond I is pro- 
portional to the modulus |cL4i| while the color represents the 
complex argument of dAi/Ai. Red: argument is +tt/2, blue is 
— 7r/2. Here dA is a gauge excitation, it corresponds to a long 
wavelength modulation of the flux through each plaquette. 
Here the wave-vector is parallel to the horizontal bonds. This 
gauge mode is associated to the smallest non-zero wave vector 
on the square lattice and is therefore four-fold degenerate. 



sq36 E=-30.93605205126 25=1 



1 2 3 4 5 



FIG. 3: Hessian eigenvector dAi corresponding to the small- 
est eigenvalue (2.62955994344) for a 36-site square lattice. 
Same representation as in Fig. 2. The complex argument of 
dAi/Ai takes only two values: ±7r/2 (blue and red), indicat- 
ing that dA is a gauge excitation. This mode corresponds to 
an increase of the flux through large horizontal loops, and no 
change for the flux going through local plaquettes. 



decreases quickly with distance (see caption of Fig. 4), 
but the fact that some non-trivial fluxes (different from 
or 7t) around the center of the excitation shows that 
it induces some non-coplanar spin-spin correlations. Far 
from the center of the excitation the spins remain in a 
collinear and in an ordered structure. The center spin 
is correlated with its neighbors, but it is very weakly 
correlated ((5*14 • Si) ~ 0) with the site which are far 
apart. Interestingly, a similar excitation is found in the 
ordered phase of the triangular lattice. 
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TABLE I: Low-energy saddle points for a square lattice cluster with 36 sites and S = 0.5. From the left column to the right: 
energy, spinon gap, lowest non-zero Hessian eigenvalue, degeneracy, number of different chemical potentials A (1 means uniform, 
etc;), minimum and maximum values of A, number of different values of \A\ minimum and maximum values of \A\, real (R) or 
complex/chiral solution (C). The spatial modulations of \Aij\ and Ai for solution a are displayed in Fig. 4. Solution b has the 
same spatial structure (two-fold degeneracy) as the one shown in in Fig. 6b. 
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FIG. 4: First excited mean field solution on a square lat- 
tice with 36 sites and S = 0.5. The modulus of the bond 
parameter (right color scale) is indicated, as well as the 
chemical potentials (left scale). This solution shows a lo- 
calized excitation (here around site number 14). The flux 
decays rapidly with distance from the center: 0.40461337T 
on plaquette [8,9,15,14], 0.037028tt on plaquette [1,2,8,7], 
-0.005291300tt on [0,1,7,6], 0.0001131tt on [4,5,11,10], etc. 
See Fig. 5 for some the spin-spin correlations in this state. 



Excited mean field solutions for k = 0.1. 
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For k = 0.1, the ground state on the square lattice is 
still uniform and without any flux, but the spinon gap 
(Fig. 1) remains finite in the thermodynamic limit. This 
mean field state has been argued b to be unstable at fi- 
nite 7Y, due to strong gauge fluctuations. The finite N 
ground state is believed to spontaneously break some lat- 
tice symmetry to form a valence-bond crystal. In fact, 
the mean field energy landscape is very different from the 
one observed in the the magnetic phase. The first local 
minima and saddle points are listed in Tab. IV A 3. 

We observed a high density of local minima and saddle 
points, with very small Hessian eigenvalues. This high 
density of excited saddle points reflects the presence of 
some strong gauge fluctuations at (large but) finite 7Y. 
We note that several such saddle points display modu- 



FIG. 5: Spin-Spin correlation in the first excited of a 36-site 
square lattice with S = 0.5 (same state as in Fig. 4). The 
radius of the circle on site j proportional to the correlation 
| (SjoSj')] with the reference site iO (the sign is available on the 
color scale) . Top: the reference spin io = 14 is at the center of 
the localized excitation. Bottom: the reference spin io = 35 
is "far" from the center of the excitation. 



lations of the bond amplitudes \A\ which have the same 
symmetries as the valence-bond crystals previously con- 
sidered for square lattice antiferromagnets. The first ex- 
citations (line ° in Tab. IV A 3) shows a one-dimensional 
modulation, which is predicted to occur at finite N when 
k = 2 mod 4.'' And among the highly symmetric so- 
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TABLE II: Low-energy saddle points for a square lattice cluster with 36 sites. S = 0.05. Solutions and b are displayed in 
Fig. 6, and solution c is shown in Fig. 7. 



lutions we also note a plaquette valence-bond crystal 
(VBC) (E = -1.82534862050) (Fig. 6). This state does 
however not correspond to a simple VBC since it is a 
complex / chiral solution with non-trivial (different from 
or 7r) fluxes on all the square plaquettes. See also Fig. 7 
for another VBC-like state. According to the analysis of 
Read and Sachdev, non-perturbative gauge fluctuations 
(proliferation of hedgehogs point-like instantons) are re- 
sponsible for the lattice symmetry breaking at finite N, 
leading to a modulation of \A\ which is exponentially 
small in N. The mechanism is somewhat different here 
since we observe VBC-like low-energy saddle points al- 
though gauge fluctuations are completely absent (A^ is 
frozen in the SBMFT). In Sec. IV C we will also report a 
case where the ground state itself exhibits some sponta- 
neous bond amplitude modulations. 
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B. Triangular lattice 

1. Hessian of the ground state 

On the triangular lattice we find that the lowest en- 
ergy state corresponds to the (spatially uniform) solution 
studied by Sachdev (0-flux state 14 ) and leads to mag- 
netic long-range order (\^3 x V3) for 2S > 0.34, and a 
gapped and deconfined Z 2 liquid for 25 < 0.34. As usual, 
these two phases can be distinguished by the spinon gap: 
it drops to zero when increasing the system size in the 
magnetic phase and stays finite in the liquid phase (see 
Fig. 8). 

The Hessian has a large lowest eigenvalue, which indi- 
cates the stability of this mean field state with respect to 
1/N fluctuations. The evolution of this lowest eigenvalue 
is plotted in Fig. 8. It slightly decreases with the system 
size (comparing 36 and 144 sites), but is certainly finite 
in the thermodynamic limit. Contrary to the square lat- 
tice situation, one does not detect any dramatic change 
of behavior between the gapped phase and the magnetic 
one. 

In fact, since the lattice is no bipartite, the ICG of 
the uniform mean field state is discrete (Z 2 ) and we do 




FIG. 6: Two low-energy mean field solutions on a square lat- 
tice with 36 sites and S = 0.05. \A\ takes only two values and 
the bonds with the stronger \A\ are indicated by fat grey lines 
(for values, see Tab. IV A 3). All the sites are equivalent and 
the chemical potential is uniform in both cases. 



not expect any gapless modes associated to small per- 
turbations Aij = A^j + dAij = with \dA l3 \ < 1. Still, 
important Z 2 (gapped) gauge excitations are expected 
in the spin liquid phase, and they will be discussed in 
Sec. IV B 3. 

The lowest eigenvector of the Hessian is represented in 
Fig. 9 for S = 0.05. As in all the cases we looked at, it 
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FIG. 7: A low-energy mean field solution on a square lattice 
with 36 sites and S = 0.05. The modulus of the bond pa- 
rameters is indicated - the right color scale shows deviations 
from the minimum value. The left color scale indicates the 
chemical potentials (dot on each site). This solution shows 
modulations with the symmetries of a 3-cell VBC. 



correspond to a phase fluctuation of the bond variables. 
The associated flux modulations for all diamond loops are 
shown in the bottom panel of the figure (note that the 
sign of each flux is somehow arbitrary since it depends 
on the choice of an origin of the loop). 
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FIG. 8: Top panel: smallest eigenvalue of the Hessian for the 
ground state state on 36-site and 144-site triangular clusters. 
Bottom: spin gap. In the thermodynamic limit this state is 
associated to magnetic long-range order (vanishing spin gap) 
for 25 > 0.34. 9 



2. Excited mean field solutions 

In Tab. Ill we list these first mean field states obtained 
for S = 0.5 in a 36-site sample. As for the square lattice 
in its magnetic phase, we only find saddle points and no 
local minima among the first states. The first excited 
saddle point is at an energy 0.94927 above the ground 
state and this gap is very likely finite in the thermody- 
namic limit. This excited state has the spatial structure 
of a localized excitations (Fig. 10) which resembles that 
of the square lattice. Finally, we note that contrary to the 
square lattice case, the lowest energy states of Tab. Ill 
have real bond amplitudes (more precisely: can be made 
real with an appropriate gauge choice), indicating the 
coplanar nature of their spin-spin correlations. 

For small enough k = 2S, the SBMFT describes a 
gapped spin liquid of Z 2 type. 9,14 Tab. IV gives the first 
saddle points obtained at 2S = 0.1. The ground state 
is uniform and all the rhombi have a vanishing flux, as 
expected. The first excited saddle point has slight modu- 
lations (~ 6.10~ 4 ) of the bond amplitudes (Fig. 11). This 
state also has vanishing fluxes on all the rhombi, but it 
differ from the ground-state by the presence of an ad- 
ditional flux 7T along some long loop winding around the 
torus (three possible choices). These three states become 
homogeneous in the thermodynamic limit and degener- 
ate with the ground state. With the ground state they 



form the four-fold topological degeneracy of the Z 2 liquid 
on a torus. ' These states were found by the optimiza- 
tion algorithm starting from random initial conditions, 
and not "forced" by hand. We are thus confident that 
the method is able to find the low energy solutions in a 
systematic way on small clusters at least up to a few tens 
of sites and bonds. 

Above the four topological ground states we observe a 
many saddle points which do not show any simple spatial 
pattern (high number of inequivalent sites and bonds). 
These states have energies significantly below the first 
vison-pair state we have found (last line in Tab. IV and 
Fig. 12). Due to the presence of complex fluxes (not or 
7r) in these states, they may not have any interpretation 
in terms of the Z 2 gauge degrees of freedom. 



3. Visons 

A Z 2 liquid possesses non-magnetic excitations named 
visons, which correspond to 7r-flux quanta of the effec- 
tive Ising gauge theory. 1 '' 5 ' 5 In a finite system (without 
boundaries) the number of vison is necessarily even, and 
a trial vison-pair state can be constructed as follows. One 
starts from the uniform mean field ground state and one 
reverses the sign of the bond parameters Aij for all the 
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TABLE III: Energy minimum and low-energy saddle points for a triangular lattice cluster with 36 sites and S — 0.5. Notice 
that the ground state is very stable (large Hessian gap: 1.051) and that all the other saddle points are unstable (negative 
Hessian eigenvalue). The gap —25.454 + 26.404 = 0.949 is also quite large compared to the kagome case (next section). The 
ground state (E = —26.40405992819) is a solution with vanishing flux on all the diamonds, and three-sublattice long-range 
spin-spin correlations (the critical value for magnetic long-range order is 2S = 0.34. '). The bond strength of the first excited 
state are displayed in Fig. 10. 
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TABLE IV: Low-energy saddle points for a triangular lattice cluster with 36 sites and S = 0.05. The ground state and the 
first excited saddle point [E — —1.76857445623, three fold degenerate) form the four-fold topological degeneracy. The other 
saddle points listed here are chiral(complex) , except for the last line, which corresponds to a pair of visons (see Sec. IV B 3 and 
Fig. 12). 



bonds ij crossing a cut extending from the a first pla- 
quette to a second one. The gauge flux is then concen- 
trated in the immediate vicinity of these two plaquettes 
which correspond to the vison core positions. Any gauge- 
invariant operator far from the vison cores is unaffected 
by this modification. 

Due to the presence of flux ir for each loop encircling 
a vison core, such a trial state is in general not self- 
consistent ((Aij) 7^ Aij and (hi) ^ k). To obtain a 
self-consistent state (saddle point) the bond moduli as 
well as the chemical potentials should be re-adjusted in 
the vicinity of the vison cores. We obtain numerically 
such self-consistent vison-pair states by looking at the 
nearest saddle point in the vicinity of the trial state (sec- 
ond stage of the algorithm described in Sec. IIIB). These 
gives access to some energetics of the visons, and to their 
energy as a function distance in particular. 

Typical vison-pair solutions are displayed in Fig. 12 
and Fig. 13. The modulation of \A\ in the vicinity of 
each vison is clearly visible, as well as the higher (less 
negative) chemical potentials in the core regions. In this 
case (S = 0.1) the vison core radius is of the order of 
two lattice spacings. These excited states are however 
not local minima, but saddle points. In all the cases 
considered here (vison distance from 2 to 4-\/3) we find 



four negative Hessian eigenvalues, with a rather weak 
dependence on the vison separation. 

A calculation of the bond modulations in a vison mean- 
field state was recently carried out by Huh, Punk and 
Sachdev ' 1 using and effective model valid in the limit 
of large spinon gap. Their result indicates that \Aij\ 
decreases on all the bonds close to the vison core. In 
our calculation it appears that the some moduli are in- 
deed depressed, but some are also enhanced (red bonds 
in Fig. 12 and grey bonds in Fig. 13). 

Some aspects of the visons energetics are summarized 
in Fig. 14. In the laige-N framework, these energies 
should be multiplied by a factor N. It can be checked 
that their mutual interaction is very weak, since the to- 
tal energy hardly depends on the distance, as expected 
in a gapped Z2 liquid phase. These calculations also 
provide some information about the spinon-vison inter- 
actions. This question is important since vison-spinon 
bound states have fermionic mutual statistics. i,! It ap- 
pears however that the spinon gap A is slightly higher in 
presence of a vison pair than in the ground state (Ao). 
This indicates some vison-spinon repulsion and makes 
unlikely the existence of a bound-state between these two 
excitations. 

We investigated the effect of inserting two ir fluxes far 
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FIG. 11: First excitation on a triangular cluster with 36 sites 
and 5 = 0.05. This state is three-fold degenerate. With the 
ground state, it forms the four-fold topological degeneracy 
expected for a Z2 liquid on a torus. Notice (scale on the right) 
that the difference between the largest \Aij\ and the smallest 
one is only ~ 6.10 -4 and should vanish in the thermodynamic 
limit. 



FIG. 9: Hessian eigenvector dA/A corresponding to the small- 
est (degenerate) eigenvalue (0.1943) for a 144-site triangular 
lattice with S — 0.05 (uniform ground state). The complex 
argument of dA/A take only two values: ±7r/2 (blue and yel- 
low), indicating that dA is a (gapped U(l)) gauge excita- 
tion. Bottom: Infinitesimal flux variation dF associated to 
the gauge mode above. For each diamond the magnitude of 
the flux variation is indicated by the radius of the dot in its 
center (see color scale for the sign). 



tri36 E=-l. 762670949826 2S=0.1 



tri36 E = -25. 45478884811 2S = 1 




FIG. 10: First excited state on a triangular cluster with 36 
sites and S = 0.5. The bond moduli are invariant by lat- 
tice rotations about the site largest chemical potential (grey 
circle). The flux vanish on all diamonds except for the six 
diamonds which diagonal bond is marked by a red dot. The 
later have flux 7r. Notice the similarity with Fig. 4. 



apart in the magnetic phase of the model. Since the 
system is magnetically ordered, visons are no longer low- 
energy excitations. Starting from a vison-pair trial state 
we look for the nearest self-consistent mean field state. A 
typical result is shown in Fig. 15, where the initial state 
was chosen to have two localized visons at the same lo- 
cations as in Fig. 13. At the end of the numerical op- 
timization (stage 2 only), it appears the algorithm has 
converged to a (unstable) saddle points where some ad- 




tri36 E=-l. 762670949826 25=0.1 




FIG. 12: Top Modulus \ A\ and chemical potential for a saddle 
point with two visons at distance d — 5v3/3, (6x6 triangular 
lattice, S = 0.1). Bottom: the flux vanishes on all diamonds 
except those with the diagonal tagged with a grey circle. The 
visons are localized on the triangles with three such circles. 



ditional pairs diamonds with flux it are present. These 
additional fluxes form an elongated ring enclosing the two 
initial vison cores. Contrary to the vison-pair states in 
the liquid phase, the spin-spin correlations are strongly 
modified all the way inside the ring (bottom panel of 
15). Indeed, the three sublattice is destroyed, although 
the spin spin correlations remain large (black circle radii) . 
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tril44 E = -49. 173695615 2S = 0.6 
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FIG. 13: Modulus \A\ and chemical potential for a saddle 
point with two visons at distance d = 4v3, (12 x 12 triangular 
lattice, S = 0.1). A string of reversed Aij (not shown) goes 
from the first vison core to the second. 
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FIG. 14: top: Energy cost of a pair of visons as a function of 
distance in a 144-site triangular cluster. Bottom: difference 
between the spinon gap A in presence of the two visons, and 
the gap Ao in the absence of visons. This difference has been 
scaled by S to compare S — 0.05 and S = 0.15). Since A > 
Ao, there is some repulsion between a spinon and a vison pair. 




tril44 E=-49. 17369561528 2S=0.6 




FIG. 15: Highly excited solution on the triangular lattice (144 
sites) with S — 0.3 (ordered phase). Top: bond amplitudes. 
Bottom: Spin-spin correlations and fluxes. If (So ■ Si) is posi- 
tive while {Si ■ Si) and {S2 ■ Si) are negative, the site i is said 
to belong to the "0" sublattice (blue circles). Likewise, the 
site which are positively correlated with site 1 (resp. 2) are 
marked in red (resp. green). The other sites are marked in 
black and the radius of the circle is proportional to the max- 
imum of I (So ■ Si)\, I (Si ■ Si) I and \{S 2 - Si}]. The diamonds 
have vanishing flux except for those marked by a grey circle, 
which have a flux ir (see text). 



This state appears to be similar to a classical vortex/anti- 
vortex pair. 



C. kagome lattice 

This section presents some preliminary results for the 
kagome lattice antiferromagnet. 

We explored the low-energy saddle points for kagome 
clusters with 36 and 48 sites, with boson densities S = 0.5 
and S = 0.1. In all these cases, the ground state is found 
to have the same modulated spatial structure. As show in 
Fig. 16, the solution spontaneously breaks some lattice 
symmetries, with two different values for the bond am- 
plitudes. This state has some magnetic long-range order 
for S = 0.5 (in the thermodynamic limit) and the three- 
sublattice structure is shown in Fig. 16. It appears that 
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FIG. 16: Lowest-energy mean field solution on a kagome clus- 
ter with 36 sites and S = 0.5. Top panel: the modulus of the 
bond parameter and the value of the chemical potentials are 
indicated. Bottom: three-sublattice structure of the spin-spin 
correlation functions. This state has correlations which are 
neither of the q = nor y3 x \/3-type. 



this state is not of the q — type, nor y/3 x Instead, 
the three-sublattice structure has a preferred lattice di- 
rection (that of lines with green and blue sites only in 
Fig. 16) and this leads to an enhancement of \Aij\ for 
the bonds aligned with that direction. The 120 degrees 
sublattice structure also implies values or tt for all the 
fluxes. In the present case the flux vanishes on all the 
hexagons, stars, and length-8 rhombi. It however has 
flux 7r on some other loops of length 8 and 10. Note 
that, as remarked by Tchernyshyov al. J> , it is not possi- 
ble to construct a state with zero flux on all loops on the 
kagome lattice. 

Our preliminary calculations on a 108 site cluster 
shows that the striped phase is no longer the ground state 
on larger systems. Due to the complex energy landscape 
(high density of low energy saddle point) and due to the 
large system size, it is numerically demanding to sample 
the solutions starting from random initial conditions on 
the 108 site kagome cluster (as we did for the 36 and 48 
sites clusters). So far we resorted to look at some known 
solutions and perturbation around them. 

On the 108 site cluster, the \fi x -\/3 turns out to be 
slightly lower than the "stripe" state (Tab. V). We con- 
clude that with 36 and 48 sites, the striped state is sta- 
bilized through resonance loops (of length 6 or 8) which 



wind around the periodic boundary conditions. Still, the 
uniform \/3 x \/3 state is not ground state either on the 
108 site cluster since since it possesses some negative Hes- 
sian eigenvalue. Work is in progress to find the absolute 
energy minimum on this cluster. 



1. Excited saddle points and Hessian eigenvalues 

A striking fact concerning the mean field states on the 
kagome lattice is the high density of saddle points and 
local minima, as can be seen in Tab. V. The energy differ- 
ences and the bond and site modulations are very small. 
For this reason, a high accuracy is needed to resolve the 
different saddle points and their (tiny) spatial modula- 
tions. It is also interesting to note that the lowest Hessian 
eigenvalues are two orders of magnitude smaller than in 
the triangular lattice. Although this should be related 
to the large degeneracy in the classical limit (large S), 
we note that the lowest Hessian eigenvalues are also sur- 
prisingly very small in the spin-gapped phase at small S 
(data not shown). This very weak energy curvature im- 
plies that, unless N is very large, the bond fluctuations 
will be very important at finite N. 

Sachdev' 1 considered the first quantum corrections to 
the energy of various coplanar and non-coplanar spin con- 
figurations. He found that coplanar state have a lower 
energy, and the degeneracy among coplanar states is it- 
self lifted (by a small amount) in favor of the y/i x \/3 
structure. His calculation assumed some homogeneous 
\Aij\ and Xi, although the self consistent states associ- 
ated to some generic classical planar state will break the 
equivalence between all the bonds and all the sites. Our 
calculations do take into account the nonequivalence be- 
tween sites and bonds and thus gives the true saddle point 
energies. The saddle points we obtain are clearly related 
classical planar states, and the underlaying planar order 
is easily revealed by looking at the spin-spin correlations. 

As an example, Fig. 17 shows the bond modulations 
and three-sublattice structure of the first excited saddle 
point (36 sites, S = .0.5). In this particular example, the 
hexagons with only two colors correspond to (slightly) 
weaker \Aij\, and thus higher energy. This is somewhat 
counter intuitive since one could have expected these 
spins to resonate between to configurations and thus to 
lower the energy on these bonds. Also, the state with 
a/3 x order is also among the low-energy local energy 
minima on the 36 site cluster (line a in Tab. V). The 
state with q = order is significantly higher and thus 
not among the first saddle points. 

Finally we note that the numerous complex (chiral) 
states we obtain must correspond to some non-planar 
spin structures. Among the non-planar states, the 
cubocl (second line in Tab. VI) is of particular interest. 
This chiral state (complex Aij) with a twelve sublattice 
magnetic structure ' 5 has recently been shown by Mes- 
sio et al. 16 to be energetically lower than the -\/3 x \/3 
state if both Aij and Bij parameters are included in the 
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FIG. 17: Excited solution on a kagome cluster with 36 sites. 
Top: bond moduli (right scale) and chemical potentials (right 
scale). Bottom: Three-sublattice structure. The radius of 
each circle is proportional to \{S a • Si)\ with a = 0, 1 or 2 
depending on the sublattice (color). 



SBMFT. The present Hessian calculation shows that it 
is a stable local minimum on the 108-site cluster, at least 
for the spin values we looked at (from S = 0.05 to 0.6), 
both in the spin-gapped (chiral spin liquid) and magnet- 
ically ordered phases. As for all the low-energy saddle 
points we found on the kagome lattice, the lowest Hessian 
eigenvalue is very small (0.0012591460 for S — 0.5 and 
1.6079020310" 5 for S = 0.05), indicating a very small 
energy scale and important finite- N corrections. 

V. CONCLUSION 

We have developed a numerical method to explore in 
a systematic way the low energy SBMFT solutions on 
finite clusters up to one hundred sites. This gives us 



access to an energy landscape which contains some infor- 
mations about excitations such has gauge modes ("pho- 
ton" in the Coulomb phase of the square lattice) and 
topological excitations such as visons (Za liquid on the 
triangular lattice). On the kagome lattice this approach 
reveals a complex landscape with tiny energy scales, and 
some unexpected magnetic ordering on small clusters. 
Work is in progress toward a better understanding of 

g kag36 E— 22,07185115576 2S-1 
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FIG. 18: Uniform mean field solution on a kagome cluster 
with 36 sites and S = 0.5. Top panel: the modulus of the 
bond parameter and the value of the chemical potentials are 
indicated. Bottom: three-sublattice structure (v3 x v3) of 
the spin-spin correlation functions. 



the large- TV kagome physics as well as quantitative evalu- 
ation of the 1 /N effects in the spin-gapped phase as well 
as in the ordered phase. 
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E 

-22.07337279228 
-22.07312100213 
-22.07185115575" 
-22.06817893625 
-22.06448876483 
-22.06414194485 
-22.06206958387 
-22.06206940890 
-22.06186658368 
-22.06162809838 
-22.06162809838 
-22.06132041734 



A 

0.12515086051 
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0.12503501657 
0.12492274319 
0.12496054795 
0.12456300988 
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0.12488134890 
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0.12509708999 
0.12484489112 



H 

0.0089249516 
0.0107154849 
0.0089311159 
-0.0112835686 
-0.0129952310 
-0.0027056846 
-0.0153407835 
-0.0161155470 
-0.0031621455 
0.0001939410 
0.0001939410 
-0.0102590748 



™°A 

-1.83200573771 
-1.83184586807 
-1.83167736298 
-1.83240200322 
-1.83195203206 
-1.83279521914 
-1.83291554283 
-1.83291658623 
-1.83448571023 
-1.83245235549 
-1.83245235550 
-1.83551900593 



max^ 
1.83113623765 
1.83157327614 
1.83167736298 
1.83097449336 
1.83036801265 
1.83033540812 
■1.83063480061 
1.83070454575 
1.83020247256 
■1.82999027588 
1.82999027588 
■1.82898234568 



2 

3 

1 

7 

7 
36 
18 

9 
18 
38 
38 
36 



0.52749717421 
0.52745562491 
0.52751942536 
0.52698803771 
0.52723842679 
0.52699979179 
0.52715425178 
0.52719208674 
0.52671374288 
0.52684964508 
0.52684964508 
0.52632645045 



m"i|A| 

0.52759396847 R. 

0.52757074339 R 

0.52751942536 R 

0.52765000379 C 

0.52764855532 R 

0.52787033175 C 

0.52763583076 C 

0.52761532928 R 

0.52785086628 C 

0.52838268073 C 

0.52838268073 C 

0.52871495375 C 



TABLE V: Some SBMFT saddle points for a kagome lattice cluster with 36 sites with 2S = K = 1. Notice also the high 
density of local minima and saddle points. The ground state (E — —22.07337279228) corresponds to Fig. 16. Line a (E = 
— 22.071851155756): uniform solution shown in Fig. 18. It has zero flux on the hexagons and 8-sites rhombi, and \/3 x \/3 
magnetic order. 
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N\ minx max^ 


N\A\ miri| A | min^i 




0-hex,0-rhom (\/3 X \/3) 
-65.19595579390" 


0.05035952577 


-0.0006338181 


1 


1 -1.82031174634 -1.82031174634 


1 0.52527760692 0.52527760692 


R 


cubocl 
-65.18108305169 


0.05032119115 


0.0012591460 


1 


1 -1.82017336946 -1.82017336946 


1 0.52524483499 0.52524483499 


C 


7r-hcx,0-rhom {q — 0) 
-64.96199404178 


0.04930776212 


-0.0660018532 


1 


1 -1.81849693150 -1.81849693150 


1 0.52476183760 0.52476183760 


R 


Stripe state 
-63.10095974957 


0.07743958702 


-31.5397600273 


3 


2 -1.78528910388 -1.72569936368 


2 0.51856978854 0.52475885718 


R 



TABLE VI: Four saddle points for a kagome lattice cluster with 108 sites with 2S = K = 1. a is uniform and corresponds 
to the Ohex-Orhom 14 state (Qi = —Q2 in Sachdev's terminology) and \/3 x \/3 magnetic r. It possesses an instability mode 
(negative Hessian value), which is a pure gauge excitation (variation of the phases of Aij only, no variation of the moduli). This 
mode is 12-fold degenerate. The cubocl was discussed using SBMFT by Messio et al. , it has a twelve-sublattice non-planar 
magnetic structure for large enough S, and it describes a gapped chiral spin liquid at small 5*. The Hessian is gapped, and the 
associated mode is a 12-fold degenerate gauge mode for 5* = 0.5. The stripe state (which is the ground state on 36 and 48-sites 
clusters) is no longer the lowest energy solution for this 108-site cluster. 
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